Finite-Size Scaling Exponents of the Lipkin-Meshkov-Glick Model 
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We study the ground state properties of the critical Lipkin-Meshkov-Glick model. Using the 
Holstein-Primakoff boson representation, and the continuous unitary transformation technique, we 
compute explicitly the finite-size scaling exponents for the energy gap, the ground state energy, 
the magnetization, and the spin-spin correlation functions. Finally, we discuss the behavior of the 
two-spin entanglement in the vicinity of the phase transition. 
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Although Lipkin, Meshkov and Glick (LMG) intro- 
duced the model baring their name in nuclear physics 
P|, it is of much broader interest. It has thus been pe- 
riodically revisited in different fields, such as statistical 
physics of spin systems 0, Q or Bose-Einstein conden- 
sates 4] to cite only a few. More recently it has also 
drawn much attention in the quantum information frame- 
work, where it has been shown to display interesting en- 
tanglement properties 0, 0, Q, Q different from those 
observed in one-dimensional models 0, 03. After al- 
most four decades, it has be en p roved to be integrable 
using algebraic Bethe ansatz [ill [l2j or mapping it onto 
Richardson-Gaudin Hamiltonians for which exact solu- 
tions have been proposed 0. However, although this 
integrability provides some important insights about the 
structure of the spectrum, it is useless to compute some 
physical quantities, such as the correlation functions, for 
a large number of degrees of freedom. 

In the spin language that we shall adopt here, the LMG 
model describes mutually interacting spins half, embed- 
ded in a magnetic field. In the thermodynamical limit, 
it undergoes a quantum phase transition that is well de- 
scribed by a mean-field analysis. This transition can be 
first or second order depending whether the interaction 
is anti-ferromagnetic or ferromagnetic. In the latter case 
and at finite size, some nontrivial scaling behavior of ob- 
servables have been found numerically 0, Q- For in- 
stance, the energy gap seems to behave as iV -1 / 3 at the 
critical point where N is the number of spins. 

In this letter, we explicitly compute these finite-size 
scaling exponents combining a l/N expansion in the 
standard Holstein-Primakoff transformation, the contin- 
uous unitary transformations, and a scaling argument. 
First, we calculate the energy gap for which we detail 
the procedure. We also give the leading finite N correc- 
tions for the ground state energy, the magnetization, and 
the spin-spin correlation functions. In a second step, we 
discuss the two-spin entanglement properties through the 
concurrence 01 which is directly related to these func- 
tions. These latter results are in excellent agreement with 
recent numerical studies 0, 0] which predict a cusp-like 
behavior of the concurrence at the transition point. 



We consider the following Hamiltonian introduced by 
LMG Q 
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where the cr Q 's are the Pauli matrices, S a = J2 i crl t /2, 
and S± = S x ± iS y . The l/N pref actor ensures that 
the free energy per spin is finite in the thermodynamical 
limit. Here we focus on the ferromagnetic case (A > 0) 
and all our results are valid for I7I < 1 (the case I7I > 
1 being trivially obtained by a simple rescaling of A). 
In this situation, the Hamiltonian @ displays a second- 
order quantum phase transition at A = \h\ 0.0] In the 
sequel, we restrict our discussion to the phase \h\ > A, 
and without loss of generality, we set h = 1. 

The Hamiltonian H preserves the magnitude of the 
total spin and does not couple states having a different 
parity of the number of spins pointing in the magnetic 
field direction (spin-flip symmetry), namely 



[H, S 2 ] = 0, and 



0, 



(4) 



for all values of the anisotropy parameter 7. In the 
isotropic case 7=1, one further has [H,S Z ] = 0, so 
that H is diagonal in the eigenbasis of S 2 and S z . Due 
to the ferromagnetic interaction between the spins, the 
ground state and the first excited state always lie in the 
subspace of maximum spin S = N/2. 

In order to analyze the spectrum of H in the large N 
limit and to capture the finite-size corrections, we per- 
form a l/N expansion of the low-energy spectrum, fol- 
lowing the ideas of Stein 0] . We first use the Holstein- 
Primakoff boson representation of the spin operator |J3 
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in the S — N/2 subspace given by 

S z = S -a f a = N/2-a^a, (5) 

S+ = (2S-a^a) 1/2 a = N 1 / 2 (l-a^a/N) 1/2 a, (6) 

S- = a) {2S - cJa) 1 ' 2 = N l / 2 a) (l~ cJa/N) 1 ' 2 ,{7) 

where the standard bosonic creation and annihilation op- 
erators satisfy [a, a'] = 1. This representation is well 
adapted to the computation of the low-energy physics 
with (a^a)/N <C 1. After inserting these latter expres- 
sions of the spin operators in Eq. ©, one expands the 
square roots in their Taylor series, and writes the result 
in normal ordered form with respect to the zero boson 
state. The Hamiltonian then reads H = Hq + + H, 
with 
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(H2) and A a — a . The index 
a keeps track of the number of bosonic operators, and 
for a given a, the superscript S codes the successive 1/N 
corrections. For instance, the nonvanishing coefficients of 



and with H 2 = 



Hq are given by h 
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Next, the Hamiltonian is diagonalized order by or- 
der in 1/N using the continuous unitary transforma- 
tion method, introduced by Wegner and indepen- 
dently by Glazek and Wilson [l9|, |20| . For apedagogi- 
cal introduction to this technique, see Ref. 21] . Note 
that the method has been applied to the LMG model in 
[22I Esl EiL I23, but its simultaneous use with the 1/N 
expansion originates in |l6j . The main idea is to diago- 
nalize the Hamiltonian in a continuous way starting from 
the original (bare) Hamiltonian H — H(l = 0). A flowing 
Hamiltonian is then defined by 



H(l) = U^(l)HU(l), 



(9) 



where I is a scaling parameter such that H(l ~ 00) is 
diagonal. A derivation of Eq. @ with respect to / yields 
the flow equation 

diH(l) = [t)(1),H(J)], where rj(l) = -U\l)diU{l). (10) 



here. This generator was first proposed in |22j, [26j and 
given a deeper physical meaning in |27j . 

More generally, to compute the expectation value of 
any operator fl on an eigenstate of H with eigen- 
value E, one must follow the flow of the operator O (Z ) = 
W{l)QU(l), by solving diQ(l) = [r](l),Q(l)}. Indeed, one 
has: 

(^\Sl\i/>) = (tp\U(l = 00) Q(l = 00) U\l = oo)|V>), (11) 

where U<(1 = oo)\ip) is simply the eigenstate of the diag- 
onal Hamiltonian H (I = 00) with eigenenergy E. In prin- 
ciple, one should follow the evolution of the S x , S y and 
S z observables, from which all others can be deduced. 
However, since we aim at computing the ground state 
magnetization and spin-spin correlation functions, and 
because of the symmetries of the model, the calculation 
can be performed more simply as follows. First, the spin- 
flip symmetry Q) implies 

(S X ) = (Sy) = 0, (12) 

(S X S Z ) = (S Z S X ) = (S v S z ) = (S z S y ) = 0. (13) 

Furthermore, since the maximum spin representation is 
one-dimensional, the coefficients of the eigenstates in this 
sector can be chosen to be real so that {{S x ,S y }) = 0. 
We are thus led to consider only the following (exten- 
sive) observables : 2S Z , 1S 2 /N, A-S 2 /N and <1S 2 /N. Of 
course, the structure of the flowing observables does not 
remain as simple as those of the initial conditions of the 
flow, even with our choice of the generator. In a nota- 
tion similar to ©, all these observables can be written 
as Q = Oo + X)fc(^fe + ^fe )' wnere the sum runs over all 
non-negative even integers fc's, and 



n = 



L^i _/ya+6- 
a,SeTN 



T and n t = J2 



( 5 ) t fc A 
jya+6+k/2-l 



(14) 

We have omitted the dependence on the flow parameter I 
of the o/s which is implicit in the following. For instance, 
the initial conditions for 2S Z are lu^ — 1, u^l = —2, 
with all other coefficients vanishing. 

The commutators [77, H] and [77, fl] are computed using 
[a, a^] = 1 and basic counting results, yielding the flows 



The anti-hermitian generator 77 (Z) must be chosen to 
bring the final Hamiltonian to a diagonal form. Wegner 
proposed r)(l) = [H d (l),H od (l)} = [H d (l),H(l)}, where 
Hd and H 0< \ are the diagonal and off-diagonal parts 
of the Hamiltonian. For our problem, it would read 
■q(l) = [H (l) , H2 (I) + H2 (l))- Such a choice suffers from 
the drawback that the tridiagonality of H(l = 0) is lost 
during the flow and that H(l) contains some terms which 
create any even number of excitations. This problem can 
be circumvented using the so-called quasi-particle con- 
serving generator rj(l) = H^il) — H^{1) that we shall use 



B,h {S) -2 V A°' n ,(S') h (8-8'+i-n) 

°l n 0,a — Z 2.^ ^o' ,a-a' -2+n n 2,a' n 2,a-a' -2+n' \ L0 ) 



n,a' ,5' 



R 0,n ,(«') h {S-S'+l-n) 

/ , u a',a-a'+n ,l 2,a l "'0,a-a'+n ' 



(16) 



n, a' ,8' 



hi s :u s : s '^l (17) 



/ j ^ l a / ,o-Q'-2+n' i '2,a' u ^2,a-a / -2+ 
n,a' ,5' 



f) X s ) _ f A k > n (S-S'+l-n) 

° lUJ k,a ~ 2^ 2, a' L' /l Q',Q-a / -2+n a; /c-i-2,Q-a / -2+n 



n,a' ,8' 



1 yyk — 2,n 
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wfc'lt^ll, (18) 
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with the definitions 

= n\(c2,C2„ 



Ak,n 
Ax' a'' 



Ca'+2C2"+/E+2 J , (19) 



a' a'' 



(20) 



C™ being the binomial coefficient a\/[n\{a — n)!]. The 
sums in 1(1511811 are constrained by the fact that all sub- 
scripts and superscripts have to be positive. For example, 
in l|15(l . n runs from to 1 + 5, a' from to a — 2 + n 
and 5' from to 5 + 1 — n. At lowest nontrivial order in 
l/N, Eqs. (PITH|l become 



0/45 - 



= -4/i 



(or 

2.0 

(0) 2 



d[h. 



2.0 



-2h (0) h {0) 



(21) 
(22) 
(23) 



These equations are the well-known Bogoliubov trans- 
form, written in a differential form [2lJ. 

In |16( , Stein integrated analytically the flow equations 
(|15fl and i|16|) for 7 = — 1, and computed the 1/A correc- 
tions to the ground state energy and to the gap. For our 
purpose, we have generalized his solution to any value of 
7, and went to order l/N 3 . Furthermore, we also com- 
puted the exact solutions to (|17fl and (|18fl and computed 
the 1/A 2 corrections to the extensive observables. De- 
tailed calculations will be presented elsewhere [2^] . 

To keep the presentation short, we will here only deal 
with the results obtained for the energy gap A. We found 



= 1 



A(A) 
A(oo) 



1 

A 



Px(A, 7 ) , Qi(A, 7 ) 



S(A, 7 )3/2 3(A, 7 ) 
P 2 (A, 7 ) , Qa(A,7) 



S(A, 7 ) 3 

f 3 (A, 7 ) 

S(A, 7 ) 9 / 2 ^ S(A, 7 ) 4 



3(A, 7) 5 /2 
, Q 3 (A, 7 ) 



(24) 



° { Nl 



where 3(A, 7 ) = (1 - A](l - 7A), A(oo) = 2 -(A^) 1 / 2 
is the mean- field gap |2j, |3j, and the JVs and Qi's are 
polynomial functions of A and 7. For the isotropic case 
7 = 1, one has Pi (A, 1) = 4A(1-A) 2 , Q X (A, 1) = -2A(1- 
A) and all contributions of order higher than l/N vanish, 
so that we recover the exact result A 7=1 (TV) = 2(1 — A) + 
2A/A. 

Let us now discuss the case 7 < 1 for which we have 
checked that A = 1 is neither a root of the Pi's nor of 
the Qi's. The result 1)24(1 shows that all l/N 1 corrections 
diverge when A approaches the critical value 1 in the infi- 
nite system, such that the larger values of i, the stronger 
the divergence. However, physical quantities cannot dis- 
play any singularity at finite TV. Using the usual ideas of 
finite-size scaling [2^ generalized in to infinitely co- 
ordinated systems, we can thus compute the scaling crit- 
ical exponents. To this end, let us suppose A close to its 
critical value 1. One can then neglect all Q terms which 



are less divergent than the P ones. In this limit, expres- 
sion 124(1 becomes a function of the variable AS(7, A) 3 / 2 , 
namely 

A(A) ~ A{oo)T A A5( 7 ,A) 3/2 , 7 , for A ~ 1. (25) 

Thus, the scaling function Pa for the gap must behave 

as [iV3(7, A) 3 / 2 ] in the vicinity of the critical point 
A = 1, for its product with the mean-field gap to be non- 
singular. Consequently, one gets A (AT) ~ A -1 / 3 . Of 
course we have only checked the scaling hypothesis up to 
the order l/N 3 , but the integrability of the LMG model 
leads us to conjecture that the very simple structure of 
the l/N expansion exhibited in ((24(1 is the same at all 
orders. 

We have performed the same analysis for the ground 
state energy, the magnetization and the two-spin correla- 
tion functions for the ground state. All results are sum- 
marized below and detailed calculations will be presented 
in a forthcoming publication |28j . 



A(A0 
eo(A) 
2(S Z )/N 
4(5 2 )/A 2 
4(5 2 )/A 2 
4<S 2 )/W 2 



aAiW 1 / 3 

-l-(l- 7 )/(2A)+a e A- 4 / 3 
1 + l/A + a 2 A~ 2/3 
a xx N- 2 ' 3 
ayyN-^ 3 

1 + 2/N + a zz N- 2/3 



(26) 
(27) 
(28) 
(29) 
(30) 
(31) 



where eo denotes the ground state energy per spin. In 
each of the above expressions, we have first written the 
(exact) non-singular contributions and, second, the term 
coming from the resummation of the most singular terms 
in the l/N expansion. Let us note however that in 1(28(1 
and 1(311) . the A -2 / 3 terms dominate the large A beha- 
vior. The coefficients a's are real numbers that cannot 
be computed within our approach since the scaling argu- 
ment only provides the exponents. Nevertheless, let us 
note that a zz — —a xx since for all A, one has 



^(<^> + <^> + <^ 2 >) = ^ = l + ^ (32) 



4S 2 

A 2 



2 

AT 



As (S* 2 ) is positive, one must also have a xx > 0. One 
can furthermore infer that N~ 4 ^ 3 corrections must exist 
in 4 (S 2 ) /A 2 and/or in 4 (S 2 Z ) /A 2 , to cancel the one of 
4(S' 2 ) /A 2 in Eq. (Eg. 

These results are in excellent agreement with the nu- 
merical data where the exponents for A and for 
4 (S" 2 ) /A 2 were conjectured. However, the scaling ex- 
ponent 2/3 for 2{S Z )/N differs from that found in [lj| 
(0.55 ±0.01). This discrep ancy comes from the too small 
system size investigated in |l5j (A = 500 spins). We have 
indeed performed a numerical study up to A = 2 14 spins 
and checked that the large A leading exponent is indeed 
2/3 0. 
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The finite-size scaling of the correlation functions also 
allows us to discuss the entanglement properties of the 
critical LMG model which have been the subject of sev- 
eral studies [ESQ Fot the ferromagnetic case consid- 
ered here, it has been shown numerically that the two- 
spin entanglement, as measured by the (rescaled) concur- 
rence [T^J, displays a singularity at A = 1. Actually, as 
shown by Wang and M0lmer [3(1] . the concurrence C for 
symmetric spin systems can be simply expressed in terms 
of the spin-spin correlation functions. More precisely, for 
the present case, one has 

(N-l)C=^{\ {Si - Si) | - N 2 /4 + (S 2 Z ) ) . (33) 



At the critical point, using the results (|29I3U|) and 
a-xx > 0, one can deduce that (S 2 — S 2 ) is positive. Then 
using l(3^|l and (|3*0jl one gets: 



(N - l)C X -_ 



N 



jN- 1/3 . (34) 



This behavior is in agreement with the numerical study 
of the finite-size scaling presented in 0, . In the ther- 
modynamical limit (N — > oo), and in the phase A < 1, 
the Bogoliubov transform i|2H23[l also gives |28| 



lim (N - 1)C A<1 = 1 - 

JV^oo \/ 1 



1-A 



7 A' 



(35) 



which generalizes to any anisotropy parameter the ex- 
pression recently given by Reslen et al. for 7 = 0. 

In summary, we have used the Holstein-Primakoff 
transformation and the continuous unitary transforma- 
tions to analyze the finite size corrections of several ob- 
servables in the LMG model. Using a 1/N expansion 
and simple scaling arguments, we have captured non- 
trivial exponents that had been conjectured since several 
decades (see e. g. 0,3) but had never found any analyti- 
cal support. This powerful combination of both methods 
clearly opens many routes to investigate. In principle, 
the physics of the broken phase (A > 1) could also be 
tackled using the same approach, after performing a ro- 
tation bringing the z-axis along one of the two directions 
of the classical magnetization. However, in this phase, 
the gap, for instance, is known to behave like exp(— aN) 
[3l| and may not be extracted from a 1/N expansion. A 
contrario, the behavior of the other quantities discussed 
here should be computed along the same line. 

Finally, we wish to underline that the results presented 
here are also relevant for the Dicke model Indeed, 
in the zero temperature limit, the LMG model can 
be put in a one-to-one correspondence with this latter 
model as recently shown in 



■cspc 
.15]. 
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